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Me»Q<i And Annnwfiis For ExtracHng Thir d Ventriri«> lrrfnr»| »ti»» 
Field of the Inven^inn 

The present invention is directed to a method and apparatus for extracUng 
third ventricle information of a brain from images thereof. 



Background of the lrn/ij>rftif»n i 
Magnetic Resonance imaging (IVIRI) can be used in diagnosis of Jarioul 
diseases in humans. The most important property to be considered in MR\ is 
the stimulation of the tissue with various radio-frequency (Rl=) pulses at 
definite time intervals and then to detect the resultant echoes. The precise 
timing of the RF pulses is vitally important for good imaging. The RF pulses 
can be repeated at a certain rate (TR) and the echoes can be detected at a 
certain time fJB), The relative time lengths of TR and TE determine the pulse 
sequences and hence the tissue visualization. 

The spin echo pulse sequence is the most commonly used pulse sequence. 
The pulse sequence timing can be adjusted to give Tl^eighted. Proton or 
spin density, and T2-weighted Images. The two variables of Interest in spin 
echo sequences are the TR and TE. All spin echo sequences include a slice 
selective 90 degree pulse followed by one or more 180 degree refbcusing 
pulses. 

A short TR and short TE will give a T1-weighted image, a long TR and short 
TE will give a proton density image, and a long TR and long TE will give a T2- 
weighted image. 

Fluid attenuated inversion recovery (FU\IR) is a type of invereion recovery 
sequence to give heavy Tl-weighting. The basic part of an inversion recovery 
sequence is a 180 degree RF pulse that inverts the magnetization followed by 



a 90 degree RF pulse that brings the residual longitudinal magnetization into 
the x-y or transverse plane where it can be detected by an RF coil. The time 
between the initial 180 degree pulse and the 90 degree pulse is the Inversion 
time (Tl). 

The spoiled gradient echo recovery (SPGR) sequence has the esame TE and 
TR as T1-werghted sequence but has an additional variable flip/Up angle of 
the spins. The flip angle is usually at or close to 90 degrees for a spin echo 
sequence but commonly varies over a range of about 10 to 80 degrees with 
gradient echo sequences. The large^ tip angles give more Tl weighting to the 
image and the smaller tip angle give mor© T2 or actually T2* weighting to the 
images. 



The size and morphology of the third ventricle is important in clinical 
pathology. As the third ventricle is situated in a very critical part deep Inside 
the brain, any lesion in the surrounding tissues would affect Its shape and 
orientation. Mass lesion in the brain would cause mass effect and directly 
influence the orientation of the third ventricle. 

Early intracerebral haemorrtiage is difflcult to visualise on CT images. The 
orientation of the third ventricle is l<ey In its identification. As tiiere is mass 
effect on one side, the third ventricle would shift from its midline position and 
its long axis would also change witii respect to the symmetry plane of the 
skull. An efficient way to extract the third ventricle plane would facilitate the 
identification of tiie early Intracerebral haemorrhage and localisation of ttie 
two landmarics. namely the anterior commissure AC and posterior 
commissure PC, for spatial normalisation of the human brain. 
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The size and width of the third ventricle are also important cHnlcal parameters. 
The third ventricle may be enlarged In. either generalised or localised 
hydrocephalus. The usual cause is blockage of the aqueduct of Sylvius^ 
Patients with Alzheimer's di8ease^ bipolar dlsoitlers' and manic depression^ 
have wider third ventricles. The width of the thW ventricle better reflects the 
degree of cholinergic deficit than the severity of histopathological changes 
such as scores of plaques and tangles in the brain of a patient with Alzheimer 
Disease^. 

Existing methods for identifying the above-mentioned pathology 
conventionally use ventricle segmentation. 

US 6 434 030 describes an automated method and/or system for identifying 
suspected lesions in a brain based on the application of a segmentation 
technique to at least one of the masked images to classify the varying pixel 
intensities and differentiate hyper-intense regions. 

US 6 205 235 illustrates a method for non-Invasive imaging of an anatomic 
tissue structure in isolation from sunounding tissues based on live-wire 
segmentation and boundary definition. 

US 6 208 347 describes a semi-automated method of MRI analysis based on 
mathematical modelling of MR\ pixel intensity histograms. 

WO g4/14132 describes a non-invasive scanning medical apparatus for 
generating an Image of at least an interior region of a subject to be examined 
The correlation of previous data to tiie scanned image is detemiined. 

IVIethods which utilise segmentation techniques can run into problems and/or 
tail when there is a serious inhomogeneity and/or noise as such systems are 
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highly vulnerable to noise. Inhomogeneity and various artefects such as 
pathology (which causes the loss of anatomical information). 

The present invention aims to substantially overcome or ameliorate the 
above-mentioned problems and the measurement of the width of the third 
ventricle will facilitate the identification of pathology. 



The method according to the present Invention allows the anatomical 
Icnowledge to be implicitly incorporated in the intelligent sampling scheme. 

The method finds application in medical imaging, in particular neutoimaging 
and provides ways for quantifying anatomicai structures. Ottier areas of 
applications include neuroinfomiatics. neurosurgery, neuroradiology and brain 
research. 

Summan/ of the irwAnti"'? 

The invention is directed to a metiiod and apparatus for quantifying the ttiird 
ventricle without segmentation and specifically, the extraction of the third 
ventricular plane and calculation of the width of the third ventricle of the 
human or animal brain in neuroimages through intelligent sampling of 
anatomical structures around the ttiird ventricle. 

According to a first aspect of ttie present Invention tiiere is provided a method 
for extracting tiiird ventricle infomnation from images of a plurality of axial 
slices of a third ventricle of a brain having an anterior commissure and a 
posterior commissure, the third ventricle having a third ventricle plane and a 
widtfi. tile method comprising: 
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a. determining a third ventricle midline for each of a number of the axial 
slices; 

b. detennining the orientation of each of the midlines; 

c. generating a histogram of the orientations of the midlines; 

d. determining the peak of the histogram to provide a peak orientation; 

e. selecting the midlines having an orientation within a predetemu'ned 
angle from the peak orientatton; and 

f. calculating the third ventricle plane from the midlines having an 
orientation within the predetermined angle from the peak orientation. 

Preferably, the step of calculating the third ventricle plane comprises 
calculating the least square fit plane of the midlines having an orientation 
within the predetemiined angle from the peak orientation. 

In a prefen-ed embodiment, the step of calculating the thiixJ ventricle plane 
further comprises: 



(I) calculating the maximum distance from the least square fit plane 
to the midlines having an orientation within tiie predetermined 
angle from the peak orientation, 

(II) generating a histogram of the maximum distance of the midlines 
having an orientation within the predetemiined angle from the 
peak orientation to the least square fit plane, 

(iii) determining the peak of the histogram of the maximum distance 
of the midlines to the least square fit plane, 

(iv) selecting the midlines lying within a predetemiined distance of 
the peak, and 

(V) recalculating the least square fit plane using the selected 
midlines to generate the third ventricle plane. 
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Preferably, the method further comprises calculating the width of the thitxl 
ventricle, by for example, determining the axial slice having the anterior 
commissure and the posterior commissure, detennlning two lines parallel to 
the third ventricle plane in said detemnined slice, said two lines being 
tangential to the image of the third ventricle in said slice to indicate the 
boundary between the third ventricle and grey matter, and calculating the 
distance between the two parallel lines, said distance being representative of 
the width of the third ventricle. 

Preferably, the step of detenninlng the thirel ventricle midline for each of a 
number of the axial slice s^ comprises calculating the local symmetry Index of 
a searching line segment, the third ventricle midline being the searching line 
segment that has the minimum local symmetry index. 

The local symmetry Index l8i(x,y.s„ e) may be calculated according to the 
following: 

|is(x,y,si.e)|xlsl(x.y.s„e) = Z 2 DifG(x8, y», s,, k) 

Oc Jf )k 

S 8 

where: 

|ls(x,y,si,e)| is the length of the searching line segment. 

Is (x.y.s,5) is the searching line segment of voxel (x.y.sO with the searching 
angle 0, and (x.y.8,) the searching point. 

cos (90" + ^ ) is denoted as c90^, 



sin (90° -i- 0) is denoted as 890^, 
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fabs (g(xs+ k X CQO0, + k x s900. s.) - g(xB- kxc905>, y, - k x q9O0, s,)) is 
denoted as DIfG (Xb. y.. 8,.k). where fabs is the absolute vaiue funcOon. the 
contributfon of voxel (Xs, y,.si) to Isi (x, y, si.^) being: 

DlfG(Xs. ye. 8,. 0.5) + DlfQ(x,. ys. s,, 1.0) + DlfG(x,. y,. s,. 3.0) + DifG(Xs, y,. s., 
5.0) +DifG(x».y„ 81.7.0). 

In a preferred embodiment, the step of determining the axial slice having the 
anterior commissure and the posterior commissure comprises: 

(1) calculating the x co-ordinate of the voxel xi for all of the axial slices 
where the third ventricle is present such that this voxel's y co- 
ordinate is the mass centre of s. yc. and (x,, yc, s,) is on the third 
ventricle plane, that is x, = -(d + c s, + b y^a. where (a. b, c) is a 
unit nomnal vector and d Is a non-positive constant; 

(2) generating the searching line segment from (xi. yc. sO such that the 
line segment is on the third ventricle plane and its centre is (x, yc. 
s,); 

(3) calculating the average grey level avg, of the searching line 
segment; 

(4) comparing the average grey level avg, for different axial slices si 
and determining the axial slice having the anterior commissure and 
the posterior commissure. 

Preferably, the step of detennining the axial slice having the anterior 
commissure and the posterior commissure comprises for T1-, PD-weighted. 
FLAIR, and SPGR MR datasets. detennining the axial slice with minimum 
average grey level avg, and for T2-weighted MR datasets it preferably 
comprises detennining the axial slice with maximum average grey level avg,. 
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According to a second aspect of the Invention there is provided apparatus 
arranged to perform a method for extracting third ventricle Information from 
images defined above. 



According to a third aspect of the invention there Is provided a computer 
program product comprising computer program Instructions readable by a 
computer apparatus to cause the computer apparatus to perfomi a method 
defined above. 



Brief Deseri ption of the r > pa v>fiiTf|f 

The present invention will now be described with reference to the sole figure. 
Figure 1, which is a flow diagram illustrating the steps Involved In an algorithm 
according to an embodiment of the present invention. 

Description of Prsferretf f Smbodimen^ 

The steps constituting a preferred embodiment of the method of the present 
invention are shown in the flow diagram of Figure 1. The method of the 
present Invention, will be discussed in more detail after a brief discussion of 
these steps. 



Given the radiological Images of the brain under consideration and the 
starting and ending axial slice (So and s„) where the third ventricle Is present 
the processing steps illustrated In the flow diagram of Figure 1 are as follows: 

SleeJ. - extract the third ventricle midline segments for all of the axial slices 
in between the starting and ending axial slices So and Sn inclusive; 



S^E^ - remove outliers of the extracted midline segments; 
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SteO - calculate the third ventricle plane (PV3) from the extracted third 
ventricle midline segment iniiers; 

~ find the axial slice (APC) In between the starting and ending axial 
slices So and s„ where the anterior commissure (AC) and posterior 
commissure (PC) are present; and 

SJgtS - in the aforementioned axial slice (APC) locate the two line segments 
parallel to the third ventricle plane (PV3) and tangential to the third ventricle, 
the distance between them is token as the width of the third ventricle. 

A brain dataset or volume is represented as a stack of parallel two- 
dimensional slices. The three dimensional volume is denoted as Vol (xyz) 
with X. y and z being the co-ordinates at voxel (x.y,2). In this case. x. y and z 
are non-negative integers satisfying 0 x Xsize. 0 y Ysize , 0 z 
Zsize where the z co-ordinate is constant on the axial slices, the y co- 
ordinate is constant on the coronal slices and the x co-ordinate is constant on 
the sagittal slices. 

If the original scanning orientation is coronal or sagittal, the axial slices are 
obtained by reorienting the original volume by reordering its voxels. The 
algorithm of the present invention worics on the axial slices. The beginning 
and ending axial slices so and s„ wheie the third ventricle is present are 
predetemiined. Any axial slice in between so and 8„ is denoted as s,. where 
s. rtself represents the axial slice as well as the axial slice number. The grey 
level at voxel (x.y.s.) is denoted as g (x,y.s,). From voxel (x.y.s.) numerous 
line segments can be drawn within s,. The line segment is denoted as Is 
(x.y.s,^) taking (x.y.s,) as its centre, with the length of line segment being a 
constant L (for example. 60 mm) and the angle with respect to the y axis 
being 0. is (x.y.s,^) is called the searching line segment of voxel (x.y.s,) 
with the searching angle 0, and (x,y.s,) is called the searching point. 
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step 1 : Extract the thtrrl wp tride midllna ? «»nrr'*P*° 

A prominent feature of the third ventricle in axial slices is that the thalamus 
(grey matter, GM) and the third ventricle (cerebrospinal fluid, CSF) are 
substantially symmetrical with respect to the third ventricle midline. On axial 
slices, the length of the third ventricle may be up to 40 mm and its width may 
vary between around 3 mm to 10 mm. The centre of the thW ventricle is 
around the mass centre of the axial slice. 

To locate the third ventricle midline in an axial slice s,. the local symmetry 
Index of a searching line segment is used to capture the anatomical futures 
of the third ventricle midline segment and thus to locate the third ventricle 
midline. Due to the variations in size of third ventricles, the local symmetry 
index should sample both the gray matter (GM) and cerebrospinal fluid (CSF). 

For the searching line segment Is (x,y,s„(9). its local symmetry index Isi 
(x.y,s,,d) measures the grey level symmetry around it. For each voxel (x..y,. 
s.) on the searching line segment, five pairs of sampling points at the opposlt^ 
sides of Is (x,y,s..6>) are taken on the lines perpendicular to Is (x.y.s,.<9) and 
passing through (xs.ys.s,) with the distance to Is (x.y,s..(9) preferably being 0.5 
mm, 1 mm, 3 mm, 5 mm and 7 mm respectively. 

cos (go» + ^ ) Is denoted as c9O0 
sin (90" + 0)\s denoted as sOO^ 

fabs (g(xs + k X c90^ , y« + ic X s90^ , s.) - g(Xs - loccOOtf , y. - k x sQOe , s.)) is 
denoted as DIfG (x«, ys, S|,lc) 

The contribution of voxel (Xs, y8.si) to Isi (x. y. S|,^) is: 

DifG(x», y„ 8„ 0.5) + DlfG(x,. y.. s,, 1.0) + DHG(x„ y,. s,. 3.0) + DlfG(x,. y„ s., 
5.0) + DlfG(X8, ys. s,, 7.0) 
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where fabsQ Is the absolute value function. 

l8i(x.y.8,,e) is the average contribution of ail the voxels on ls(x,y.s,.e). that Is. 

|te(x,y,s,,e)|xl8i(x.y,s,.e)s:S 2DifG(x„ys.8,.k) 

(X .y }k' 

8 8 

Where |ls(x.y.s,.G)| is the length of the searching line segment in millimeters 
(mm). 

The third ventricle midline segment on axial slice s, Is the searching line 
segment that has the minimum local symmetry index. The extracted third 
ventricle midline segment is called the approximated third ventricle midline 
segment (ATS/IWIS). 

Step ^\ Rgpnoye outliers of thfl extraftteri miriii p a saomante 
The approximated third ventricle midline segments (ATVIWlSs) are processed 
In two steps, to remove outliers. In the manner described for example In the 
applicants copending Intematlonal Patent Application PCT/SG02/00231, the 
content of which is incorporated herein by way of reference. 

Firstly, the orientations of ail the ATVMSs are calculated and a histogram of 
the orientations is obtained. The peak of the histogram Is determined and is 
called the peak orientation. Those ATVMSs with an orientation deviating 
from the peak orientation by more than a predetermined value, for example 
r, are considered as orientation 'outiiere' while the rest of the ATVMSs ai« 
considered to be orientation 'iniiers'. 

Secondly, the least square fit plane of the orientatton iniiers is calculated. 
The maximum distance of ail the orientation iniiers to this plane Is calculated 
and the peak of the histogram of all the distances is obtained. Those 
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Orientation inliere with a distance deviating from the peak distance by more 
than a value of. for example 1mm. are considered the third ventricle plane 
outliers, while the rest of the orientation inliers are considered as the third 
ventricle iniiers. 



step 3: Catailate the third w^nfrf i^if p|f, n -> 

The third ventricle plane is approximated ftom the third ventricle iniiers using, 
for example, the least square fit plane of the third ventricle iniiers. The third 
ventricle plane is denoted as: 
ax + by + cz + d=!0 

where (a. b. c) is a unit nomial vector and d is a non-positive constant 

Step 4: Find the axial slicf» with the anlftrl«|r a n d ttosterior m,n >v|f««.fr»o 
Any method for Identification of the anterior commissure (AC) and posterior 
commissure (PC) may be used to locate the axial slice with the two 
commissures thereon (APC). This may also be identified In the following 
way: 

1 . Calculate the x co-ordinate of the voxel x, for all of the axial slices s, In 
between the beginning and ending axial slices so and s„ where the third 
ventricle is present such that this voxel and the mass centre of s. have the 
same y coordinate yc, and (x,. yc. sO is on the third ventricle plane, that is x, = - 
(d + c si + b yc)/a. 

2. Form the searching line segment from (x,. yc. s^ such that the line 
segment is on the third ventricle plane and its centre is (x, yc. s,). 

3. Calculate the average grey level of the searching line segment. For 
the axial slice s„ the calculated average grey level is denoted as avg,. 

4. Compare the average grey level avg, for different axial slices s,. i^or 
T1.. PD-welghted. FU^IR. and SPGR IVIR datasets. the axial slice with 
minimum avg, is taicen as APC. For T2.weighted MR datasets. the axial slice 
with maximum avgr is taiten as APC. 
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Calculate the third ventricle width by locating the left-most and right- 
most lines parallel to the third ventricle plane and tangential to the third 
ventricle in the APG. that Is the boundary between the third ventricle and the 
grey matter. The distance between the two parallel lines is defined as the 
third ventricle width. 



In summary, the present invention Is directed to a method of extracting the 
third ventricle plane which is robust to noise, inhomogenelty and various 
artefacts. It is also directed to calculating the width of the third ventricle of a 
brain from neuro images. 

Extracting the third ventricle plane and measuring the width of the third 
ventricle is of clinical importance for both pathology detection and 
morphological description of brains. The present Invention- proposes a fast 
and automatic method for quantifying the third ventricle based on intelligent 
sampling of anatomical structures, namely the thalamus and the third 
ventricle, around the third ventricle based on the combination of anatomical 
knowledge and image analysis technique. 

In contrast to conventional methods in which the third ventricle is segmented 
the method embodying the present invention extracts the midlines of the third 
ventricle based on the local symmetry of the cerebrospinal fluid (the third 
ventncle) and the grey matter (the thalamus). The third ventricle plane is 
taken to be the least square fit plane of ail the midlines of the third ventricle 
The width of the third ventricle is calculated as the distance between two lines 
parallel to the third ventricle plane and tangential to the third ventricle on the 
axial slice containing the anterior and posterior commissures. 
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Claims: 

1 . A method for extracting third ventricle Information from images of a 
plurality of axial slices of a third ventricle of a brain having an anterior 
commissure and a posterior commissure, the third ventricle having a 
third ventricle plane and a width, the method comprising: 

(a) determining a third ventricle midline for each of a number of the axial 
slices; 

(b) determining the orientation of each of the midHnes; 

(c) generating a histogram of the orientations of the midlines; 

(d) determining the peak of the histogram to provide a peak orientation; 

(e) selecting the midlines having an orientation within a predetemiined 
angle from the peak orientation; and 

(f) calculating the third ventricle plane from the midlines having an 
orientatton within the predetermined angle from the peak orientation. 

2. A method according to claim 1 wherein the step of calculating the third 
ventricle plane comprises cateulating the least square fit plane of the 
midlines having an orientation within the predetermined angle from the 
peak orientatton. 



A method according to claim 2 wherein the step of calculating the third 
ventricle plane further comprises: 

(I) calculating the maximum distance from the least square fit plane 
to the midlines having an orientation within the predetermined angle 
from the peak orientation. 

(ii) generating a histogram of the maximum distance of the mkilines 
having an orientation within the predetermined angle from the peak 
orientation to the least square fit plane. 
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(Hi) determining the peak of the histogram of the maximum distance 
of the midlines to the least square fit plane, 

(iv) selecting the midHnes lying within a predetennined distance of 
the peak, and 

(V) recalculating the least square fit plane using the selected 
midlines to generate the third ventricle plane. 

A method according to any one of the preceding claims, .further 
comprising calculating the width of the third ventricle. 

A method according to claim 4. wherein the step of calculating the 
width of the third ventricle comprises determining the axial slice having 
the anterior commissure and the posterior commissure, detemiining 
two lines parallel to the third ventricle plane in said determined slice, 
said two lines being tangential to the image of the third ventricle in said 
slice to indicate the boundary between the thirel ventricle and grey 
matter, and calculating the distance between the two parallel lines, said 
distance being representative of the width of the thirel ventricle. 

A method according to any one of the preceding claims, wherein the 
step of determining the third ventricle midline for each of a number of 
the axial slice s, comprises calculating the local symmetry Index of a 
searching line segment, the third ventricle midline being the searching 
line segment that has the minimum local symmetiy Index, 

7. A method according to claim 6. wherein the local symmetry index 
lsi(x,y,8|, 6) Is calculated according to the following: 



6. 



|ls(x.y,si,e)IxlsKx.y,s„9) = Z S DifG(x8, ys. s„ k) 
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where: 

|l8(x,y,si,e)| is the length of the searching line segment. 

Is (x.y,s,tf ) is the searching line segment of voxel (x,y.si) with the searching 
angle 9, and (x.y.si) the searching point, 

cos (go* + ^ ) Is denoted as c90d, 

sin (90° ■<■ ^) Is denoted as s9O0, 

febs (g(xs + k X c90^ . y, + k X 890^ . s.) - g(X8 - kxc90<9 . y, - k x s90^ . s,)) is 
denoted as DifQ (x,. ys. s„k), where f^bs Is the absolute value function, the 
contribution of voxel (x,, y„si) to Isi (x, y, s..^) being: 

DifG(x». ys. s,, k1) + DifG(Xs. ys. s,. k2) + DifG(Xs. y,. s,. k3) + DlfG(xs, ys. s,, k4) 
+ DifG(Xs, ys. 8„ k5). k1. k2. k3, k4, and kS are constants. 

8. A method according to claim 7, wherein k1 is around 0.5mm. 

9. A nr>ethod according to claim 7, wherein k2 is around 1mm. 

10. A method according to claim 7, wherein k3 Is around 3mm. 

11. A method according to dalm 7, wherein k4 is around 5mm. 

12. A method according to claim 7, wherein kS is around 7mm. 

13. A method according to dalm 5, wherein the step of determining the axial 
slice having the anterior commissure and the posterior commissure 
comprises: 

(1) calculating the x co-ordinate of the voxel x, for all of the axial 
slices where the third ventricle Is present such that this voxel's y co- 
ordinate is the mass centre of s, yc. and (x,. y,. s,) is on the third 
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ventricle plane, that Is jq « + c s, + b y^/a, where (a. b. c) Is a unit 
normal vector and d Is a non-positive constant; 

(2) generating the searching line segment Irom (x„ yc, s,) such that 
the line segment is on the third ventricle plane and Its centre Is (x, yc 

si): 

(3) calculating the average grey level avg, of the searching line 
segment; 

(4) comparing the average grey level avg, for different axial slices s, 
and detemilning the axial slice having the anterior commissure and the 
posterior commissure. 

14. A method according to claim 13 wherein the step of detemilning the axial 
slice having the anterior commissure and the posterior commissure comprises 
for T1-. PD-welghted. FU^IR. and SPGR MR datasets, detemiining the axial 
slice with minimum average grey level avgi. 

15. A method according to claim 13 wherein the step of detemiining the axial 
slice having the anterior commissure and the posterior commissure comprises 
for T2-weighted MR datasets comprises detemilning the axial slice with 
maximum average grey level avg..16. An apparatus arranged to perform a 
method for extracting third ventricle information from images according to any 
one of the preceding claims.17. A computer program product comprising 
computer program instructions readable by a computer apparatus to cause 
the computer apparatus to perfomi a method according to any one of claims 1 
to 15. 
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ABSTRACT 

■Method And Appqmtiis For Extrarfinn Thi rd Ventrirte lnfnr.y, ^*>»,^ 

A method for extracting third ventricle Information from Images of a plurality of 
axial slices of a third ventricle of a brain comprises determining a midline for 
each of a number of the axial slices, determining the orientation of each of the 
midlines, generating a histogram of the orientations of the midlines, 
detemiining the peak of the histogram to provide a peak orientation, selecting 
the midlines having an orientation within a predetermined angle from the peak 
orientation and calculating the third ventricle plane from the midlines having 
an orientation vi^lthln the predetermined angle from the peak orientation. 



Figure 1 
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1. Extract the third ventricle midlines 
"I 



2. Remove outliers 



3. Approximate the third ventricle plane 



4. Find the axial slice with AG and PC 



5. Locate the two line segments parallel to 
the third ventricle plane and tangent to 
the third ventricle 



